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ABSTRACT 

In this Letter we investigate the origin of the near-infrared emission of the Herbig Ae star MWC 758 on sub-astronomical 
unit (AU) scales using spectrally dispersed low resolution (R=35) AMBER/VLTI interferometric observations both in 
the H (1.7/im) and K (2.2/im) bands. We find that the K band visibilities and closure phases are consistent with 
the presence of a dusty disk inner rim located at the dust evaporation distance (0.4 AU) while the bulk of the H 
band emission arises within 0.1 AU from the central star. Comparing the observational results with theoretical model 
predictions, we suggest that the H band emission is dominated by an hot gaseous accretion disk. 
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1. Introduction 

The development of long baseline near-infrared (NIR) in- 
terferometry allows us to spatially resolve the emission of 
the innermost regions of disks around young pre-main se- 
quence stars (<1 AU in the nearest star forming regions). 
These regions are important for understanding the star-disk 
interaction, the gas accretion process onto the star, and the 
launching of jets and winds. Moreover, planet formation is 
likely to be favored in the high density regions of the inner 
disk. 

The existing NIR interferometric observations of low 
mass T Tauri and intermediate mass Herbig Ae disks (HAe 
! hereafter; Millan-Gabet et al. 2007, Monnier et al. 2006, 
Akeson et al. 2005, Eisner et al. 2004, 2007), show that 
emission in the K band (2.2 /im) is dominated by dust at 
the sublimation temperature, located in a puffed-up rim 
(Dullemond et al. 2001, Isella and Natta 2005). For the 
more luminous HBe stars, dust emission alone can not ex- 
plain the observed visibilities, and it has been suggested 
that the K band excess flux is dominated by the emission of 
a dense, optically thick gaseous disk which extends well in- 
side the dust sublimation radius (Malbet et al. 2007) . Hints 
of the presence of an hot gaseous disk are also provided 
by recent spectrally dispersed low resolution observations 
of a sample of HAe stars where dust emission alone can- 
not reproduce the spectral variation of the visibility across 
the K band (Eisner et al. 2007). Most gas tracers, such 
as the hydrogen recombination lines, are very likely emit- 
ted by matter which is either accreting or ejected from the 
system (Malbet et al. 2007; Tatulli et al. 2007; Kraus et 
al. 2008; Eisner et al. 2007b), rather than from the gaseous 
disk itself. In only one case (51 Oph; Tatulli et al. 2008) 
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Amber/VLTI medium spectral resolution observations of 
the CO v=2-0 lines at 2.3 fim confirm their origin in the 
inner disk, inside the dust sublimation radius. 

In this Letter we present spectro-interferometric obser- 
vations of the HAe star MWC 758 (ST A5, D=200 pc, 
L»=22 L Q , M»= 2 M Q ) across the H and the K band ob- 
tained with AMBER/VLTI. We show that the H and K 
band excess emission arises from two different disk regions, 
the first of which is much more compact and hotter than 
the second. We identify the dominant source of the 2 /im 
emission as the dust rim at the sublimation radius, and the 
dominant source of the 1.7 /im emission as the gaseous disk 
inside this radius. 

The structure of this letter is as follows. Sec. 2 describes 
the observations and the data reduction, Sec. 3 describes 
the results, which are discussed in Sec. 4. Sec. 5 summarizes 
briefly our conclusions. 

2. Observations and data reduction 

MWC758 (K — 5.8 mag, H = 6.5 mag) was observed 
with AMBER on December 30, 2006, and on January 
8, 2007, with three 8 m Unit Telescopes (UTs) of the 
VLTI, using the UT1-UT3-UT4 triplet. On both dates, the 
source was observed close to transit; the hour angles (0.25 
and 0.54 hours, respectively) of the two sets of observa- 
tions were separated by 20 min. At the time of the ob- 
servations, the parameters of the three projected baselines 
were: B=81m, PA=45° (UT1-UT3), B=59m, PA=102° 
(UT3-UT4), and B=123m, PA=69° (UT1-UT4), on Dec.30, 
2006, and B=84m, PA=45° (UT1-UT3), B=58m, PA=101° 
(UT3-UT4), and B=125m, PA=68° (UT1-UT4), where the 
position angle PA is measured east of north. Observations 
were performed in low spectral resolution mode (R=35), in 
both the H and K bands. 

Raw visibilities and closure phases as a function of wave- 
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gas emission 



T e// = 2500 K 
Ri„ > 0.05 AU 
Rout <0.1 AU 
X„ = 0.41 
X K = 0.33 



dust 



emission 



X H -- 
X K = 



~ 1400 K 
= 0.34 AU 
= 0.39 AU 
= 0.50 
= 0.79 



Rout 

x H = 



- 1300 K 
: 0.40 AU 
= 0.49 AU 

0.18 

0.46 



incl = 40° 



incl = 30° 



Table 1. Parameters of the best fit models shown in Fig.Q] 
andO Model I corresponds to the rim model described in 
Isella and Natta (2005); in Model II we have added to the 
star and rim an additional source of emission, as described 
in the text. 



length were computed followin g the standard da ta reduc- 
tion algorithms described in Tatulli et al.(2007)j for each 
spectral band separately. The raw visibilities were corrected 
for atmospheric and instrumental effects by dividing them 
by the raw visibilities measured on a reference source. As 
a calibrator, we chose HD34053 (K = 5.9, H — 5.9) which 
has similar magnitudes and, according to the CADAR cat- 
alog (Pasinetti Fracassini et al. 2001), has a diameter of 
0.2 mas - that is, unresolved by the interferometer. The 
closure phase is theoretically independent on the Earth's at- 
mosphere (see, e.g., Monnier et al. 2007 for the definition of 
closure phase) and does not need to be calibrated. However, 
AMBER introduces an offset in the closure phase which is 
corrected by subtracting the calibrator closure phase from 
the MWC758 observed values. 



3. Results 

The results are presented in Fig. [TJ which shows the vis- 
ibilities and closure phases as function of wavelength for 
the two different data sets: panel (a) refers to the Dec. 30, 
2006 observations (hour angle HA=0.25 h), panel (b) to 
the Jan.8, 2007 one (HA=0.54 h). The results from the two 
data sets practically coincide, as expected given the small 
difference in HA. Due to the lower flux in the H band, data 
between 1.6 and 1.8 /Ltm are averaged in order to increase 
the signal to noise ratio. Even so, the error on the H band 
visibility is about 30%, much larger than that of the K band 
data (~5%). 

In the K band, the source is resolved on all the base- 
lines, which provide an angular resolution of 1.8 mas, cor- 
responding to 0.35 AU at the distance of MWC 758. The 
closure phase is always larger than zero and increases with 
wavelength between 2 and 2.5 /im. Note that a non zero 
closure phase implies that the K band emission is not cen- 
trally symmetric. This measurement alone rules out models 
such as spherical envelopes, flat disks, etc. for the K band 
emission (see also the discussion in Monnier et al. 2006). 

In the H band, the MWC758 emission is unresolved 
on the UT3-UT4 and UT1-UT3 baselines, with a visibil- 
ity equal to 1 within the error bars; on the longest baseline 
UT1-UT4 the visibility is slightly smaller than unity for 
both data sets (V= 0.8 ± 0.1). The H band closure phase 
is consistent with 0, as expected given that the emission is 
practically unresolved. 



4. Discussion 

To analyze the AMBER data, we compute visibilities and 
closure phases as a function of wavelength using the rim 
model of Isella and Natta (2005), where a complete descrip- 
tion of the physical models and of their limitations can be 
found. Here we just want to remark that, once the stellar 
parameters are defined, the rim location and shape depend 
only on the grain emissivity, i.e., for known dust compo- 
sition, on the grain size distribution. We assume silicate 
grains, which are the grains with the highest sublimation 
temperature according to Pollack et al. (1994). The fitting 
procedure is described in Isella et al. (2006) and the rim 
parameters that give the best fit to the data are given in 
Table 1. Ri„ is the distance of the rim from the star on 
the disk midplane, R ou t the distance of the top of the rim, 
so that AR=R ollt — Ri„ measures the curvature of the rim 
due to the dependence of the dust sublimation temperature 
on the gas density. The rim emission can be approximately 
described as a black-body at T e ff- The two quantities X# 
and X^f give the fraction of the dusty rim and the gaseous 
disk emission (see the discussion below) to the total flux in 
the H and K band respectively; the missing flux to reach 
the unity is due to the stellar photosphere. 



4.1. Star + rim models 

Model I includes the star and the rim emission only. As 
shown in Fig. [1] (solid curves), it provides an adequate, 
although not perfect, description of the K band visibility 
and closure phase measurements. Both the derived rim in- 
ner radius (0.34 AU) and the disk inclination (40°) are in 
agreement with the values obtained in Isella et al. (2006) 
based on earlier PTI observations; the disk position angle 
is not well constrained by our data, varying between 90° 
and 150°. 

However, such models fail to account for the observa- 
tions in the H band, where they predict that the total emis- 
sion should be spatially resolved. The reason is clear if we 
consider that the rim is physically narrow, e.g. AR ^0.05 
AU, and has practically the same size at all wavelengths. 
In this Ccl.SC j i\ variation of visibility and closure phase with 
wavelength may occur for two reasons only, one intrinsic, 
because the contribution of the unresolved stellar flux in- 
creases at shorter wavelengths, and one instrumental, since 
the interferometer resolution decreases with wavelength. In 
the calculation of the visibility these two effects tend to 
compensate one other and the resulting values are roughly 
constant, contrary to what we observe. The discrepancy 
between observations and model predictions is even more 
striking if we consider the dependence of the closure phase 
on wavelength. The higher angular resolution achieved in 
H band amplifies the asymmetry of the model emission, re- 
sulting in a very large values of the closure phase in the H 
band, contrary to what is observed. Note that a discrep- 
ancy between the rim model and the observations is also 
present in the SED (Fig. [2] panel (a)) since the best fit rim 
model can account only for about the 80% of the emission 
measured in the H band. 



4.2. Evidence for gas emission in the inner disk 

Motivated by these considerations, we modified our rim 
model introducing an additional source of emission, which 
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Fig. 1. Spectrally dispersed visibility and closure phase measured across the H (1.65 /zm) and K (2.2 /mi) bands. Panels 
(a) and (b) show the data for two different hour angles, HA=0.25 h and HA=0.54 h respectively. The three upper panels 
show calibrated visibilities for the three different baselines as labeled in figure and described in Sec. 2; corresponding 
closure phases are plotted in the bottom panel. Visibility errors across K band are between 3 and 5%; H band data has 
been spectrally averaged to increase signal to noise ratio, the associated error is about 30%. Solid lines correspond to the 
rim model (I), dashed lines to the gas+rim model (II). Model parameters arc reported in Tab. 1. 



we describe as a uniform black body ring where the ring 
radii and its effective temperature are free parameters. The 
values of the parameters of the best-fitting model (Model 
II) are given in Table 1. As shown in Fig. [T] and Fig. [2] 
(dashed lines) a good fit to both the H and K band visi- 
bilities and fluxes is obtained for T e ff = 2500 K and ring 
outer radius <0.1 AU. The exact values of the ring radii are 
not well constrained since the H band emission is poorly re- 
solved. However, the lower H band visibility measured on 
the longest baseline seems to indicate ring radii between 
0.05 and 0.1 AU. Fig. 2(b) shows that this model fits very 
well the MWC 758 SED. Such a model is clearly a simplified 
representation of a more realistic inner disk in which the gas 
temperature and emission resonably decrease smoothly for 
increasing stellocentric distances. Nevertheless, it can re- 
produce the rough spatial scale for the gas emission and 
it is therefore adequate to understand the current observa- 
tions. 

The physical origin of this additional emission is hard 
to constrain from the available data. We suggest that it 
can be identified with the emission of the gaseous accre- 
tion disk inside the dusty rim, based on the comparison 
of the observational constraints we have with the predic- 
tions of the Muzerolle et al. (2004) models. These authors 
model the physical structure of the inner disk of a typical 



Herbig Ae star, including dust evaporation and the forma- 
tion of the rim. Inside the rim, the gaseous disk proper- 
ties depend on the mass accretion rate M acc . We estimate 
that in MWC 758 Af QCC - 2 x 10~ 7 M Q /yr, based on the 
strength of the Br7 emission^, as described in Garcia Lopez 
et al. (2006). For such an M acc , the models predict that the 
gas in the disk mid plane reaches a temperature of about 
2500 K at a distance of 0.1 AU, consistent with our con- 
straint that the hot component is located within 0.1 AU 
from the star. Fig. [3] shows the emission of the gas disk 
for two accretion rates (10 -7 and 10 -6 M Q /yr), which en- 
compass the estimated value in MWC 758. The 2500 K 
black body emission is roughly intermediate between the 
two, and its spectral shape is similar, supporting the inter- 
pretation of the "second ring" as gas emission within the 
dust sublimation radius. 

One should note that the Muzerolle et al. (2004) models 
have been computed for slightly different stellar parame- 
ters, so that the comparison is necessarily very qualitative. 
Also, we were not able to compute the expected visibilities 
and phase closures of the Muzerolle et al. (2004) models, 
since only the wavelength dependence of the spatially inte- 
grated flux is available. 

1 Telescopio Nazionale Galileo unpublished data 
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Fig. 2. SED of MWC 758. Symbols show the photometry from the literature (Eisner et al. 2004, Malfait et al. 1998, 
Bcskrovnaya et al. 1999, Cutri et al. 2003) corrected for an extinction Ay = 0.4, R=3.1. The stellar spectrum (short- 
dashed line) is that of a Kurucz model with T e //=8000, Log </=4. Panel (a) shows the flux predicted by model I (solid 
line) composed by the dusty disk rim emission (short-long dashed line) and the star. Panel (b) shows the prediction for 
model II where an additional black body component is added (dashed line) . In each panel the image of the corresponding 
model is shown in the inset. Refer to Tab. 1 for model parameters. 




-9.5 



0.2 
log(X) (urn) 

Fig. 3. Comparison of the 2500 K black body emission, 
added in model II (solid line), with the emission of the inner 
gaseous disk predicted by Muzerolle et al. (2004) for two 
mass accretion rates, 10 -6 M Q /yr (top curve) and 10~ 7 
Mo/yr (lower curve), over the relevant reange of wave- 
length between 1 and 2.5 /imrespectively. 



We finally remark that the introduction of the addi- 
tional black body emission has consequences for the rim 
structure and for the derived dust grain size. In particu- 
lar, going from model I to model II, the inner rim radius 
increases from 0.34 AU, obtained with a dust grain size of 
1 fim, to 0.40 AU corresponding to dust grains of 0.5 jum. 
This confirms (Isella et al. 2006) that the grains in the in- 
ner disk of MWC 758 are on average larger than the typical 
dust size, i.e. 0.01 /jm, of grains in the interstellar medium. 
However, an accurate determination of the grain properties 
requires to understand the global properties of the inner 
disk, gas and dust. 



5. Conclusions 

In this Letter we present spectro-interferometric observa- 
tions (R=35) of the pre-main sequence HAe star MWC 758, 
obtained with VLTI/ Amber in the H and K bands. The ob- 
servations have been performed with three UT telescopes 
in order to measure both spectrally dispersed vibilities and 
closure phases. The H-band data and the closure phase are 
essential quantities in order to distinguish between differ- 
ent models of the disk structure. The observations provide 
a spatial resolution of 0.8 mas in the H band and 1.8 mas 
in the K band, corresponding to scale sizes of 0.16 and 0.35 
AU, at the distance of MWC 758. 

In the K band, the source is resolved with non-zero clo- 
sure phase (about 10°), as expected if the emission is due 
to dust in a rim at the sublimation radius (0.49 AU), seen 
at an inclination of about 30°. In the H band the source is 
only barely resolved with the longest baseline (125 m), im- 
plying that the emission at these wavelengths arises within 
0.1 AU from the central star. 

We suggest that moving in wavelength from ~2.5 ^m 
to the H band, centered at 1.7 fim, one moves from rim- 
dominated emission to a hotter, more compact component 
which is qualitatively consistent with the emission of the 
dust-free inner gaseous disk, as modeled by Muzerolle et 
al. (2004). These results support the suggestion of Eisner 
et al. (2007) that an extra component in addition to the 
dust would improve the fit to the spectrally dispersed, low 
resolution K band visibilities obtained with the Palomar 
Testbed Interferometer for a small sample of HAe stars. 

Our results show that the innermost gaseous disk can 
be detected and studied with current NIR interferometers 
through the detection of the continuum gas emission in the 
H band, where the gas contribution is higher than the stel- 
lar and dust emission, which dominate at shorter and longer 
wavelengths, respectively. Acccurate models of the spatial 
and spectral dependence of the inner gaseous disk emission 
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for a range of stellar and accretion properties, are much 
needed. 
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